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ABSTRACT 

We report on the results of an i-band time-series photometric survey of NGC 2516 
using the CTIO 4m Blanco telescope and 8k Mosaic-II detector, achieving better than 
1% photometric precision per data point over 15<i<19. Candidate cluster mem- 
bers were selected from a V vs V — I colour magnitude diagram over 16 < < 26 
(covering masses from 0.7 down to below the brown dwarf limit), finding 1685 
candidates, of which we expect ~ 1000 to be real cluster members, taking into ac- 
count contamination from the field (which is most severe at the extremes of our mass 
range). Searching for periodic variations in these gave 362 detections over the mass 
range 0.15 < M/M© < 0.7. The rotation period distributions were found to show a 
remarkable morphology as a function of mass, with the fastest rotators bounded by 
P > 0.25 days, and the slowest rotators for M < 0.5 Mq bounded by a line of P oc M^, 
with those for M > 0.5 Mq following a flatter relation closer to P '--^ constant. Models of 
the rotational evolution were investigated, finding that the evolution of the fastest rota- 
tors was well-reproduced by a conventional solid body model with a mass-dependent 
saturation velocity, whereas core-envelope decoupling was needed to reproduce the 
evolution of the slowest rotators. None of our models were able to simultaneously 
reproduce the behaviour of both populations. 



Key vifords: open clusters and associations: individual: NGC 2516 
photometric - stars: rotation - surveys. 



techniques: 



1 INTRODUCTION 

NGC 2516 is a well -studied nearby Interm ediate-aKe open 
cluster [r^ 150 Myr; IJeflFries. Thurston fc H amblv 2001), at 
a dist ance modulus (M — m)o = 7.93 ±0.14 (jTerndrup et al.l 
l2002h . The canonical reddening along the line of sight to the 
cluster is E{B — V) — 0.12, implying a pho tometrically de- 
rived metallicity of [Fe/H] = -0.05 ± 0.14 (|Terndrup et all 
L2OO2). The same authors derive a spectroscopic metallicity 
from two stars of [Fe/H] = 0.01 ± 0.03, consistent within 
the errors, and implying a slightly larger distance modulus 
(M — m)o — 8.05 ± 0.11. We have adopted the smaller dis- 
tance modulus of (M — m)o = 7.93 ± 0.14 and photometric 
metallicity for the remainder of this work, but any differ- 
ences induced in the results by the uncertainty in distance 
modulus are nevertheless small. 

NGC 2516 has often been referred to as the 'South- 
ern Pleiades' owing to its similarity both in angul ar extent 
and age, where we adopt an age of 100 Myr (|Mevnet et al.l 



for the Pleiades, using the classical main sequence 
turnoff ages throughout this work for consistency with the 
estimate for NGC 2516. The Pleia des metallicity is very sim - 
ilar, at [Fe/H] = -0.034 ± 0.024 (|Boesgaard fc Frie]||l990l '). 
Clu sters in the age range between a Persei (~ 50 Myr, 
e.g. |l3 asri fc Martin' '1999') and the Hyades (~ 625 Myr, 
IPerrvman et al. 1998) arc important for constraining stellar 
evolution models. Figure [T] shows a plot of the age at which 
a solar metallicity star reaches the zero age main sequence 
(ZAMS) as a function of stellar mass, indicating that solar 
mass stars reach the ZAMS at ~ 40 Myr, and ~ 0.4-0.5 Mq 
stars reach the ZAMS at the age of NGC 2516. It is there- 
fore highly useful in this regard owing to both its age and 
relative proximity, where we can study solar mass stars after 
evolving for a short time on the main sequence, and lower 
mass stars as they finish pre main sequence evolution. The 
similarity in age to the Pleiades (modulo the uncertainty 
in this parameter) means that comparisons of these clusters 
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Figure 1. Plot of the approximate age at which the ZAMS is 
reached for solar m etallicity stars, from the NextGen models 
tearafFe et al.|[l998l) . computed as the time taken for the star 
to reach within 1% of its 1 Gyr age radius. The horizontal dashed 
line indicates the approximate 150 Myr age of NGC 2516. 



could be used to study the dependence of stellar parameters 
on environment. 

NGC 2516 is very well-studied in the literature. The 
most relevant existing work includes me mbership sur- 
veys IjJeffries et al ] I2OOII . ISung et all I2OO2I ). characterisa- 
tion of the mass function from solar mass llJeffries et al 
2001) down to the brown dwarf regime (|Moraux et al 



20051 ). a survey for radial velocities an d binarit y at the 



bright end of the cluster sequence ({Gonzalez fc Lapassetl 
I2OO0I ). and studies of rotation using vsini covering stars 
down to G spectral types |Terndrup et al.l |2002| and 
I Jeffries. James fc ThurstonlllQQa ). A number of X-ray sur- 
veys, including monitoring for variabi lity, have al so been car- 
ried out (e.g. Wolk et al. 2004, Ramsay. Harra fc Kav.20o3 . 



iDamiani et al.f 20031 iSciortino et al.ll200lh , the most recent 
of which (|Pillitteri et all 120061 ) used a deep XMM-Newton 
observation to probe the low-mass members (down to ~ M5) 
of the cluster. They found that NGC 2516 stars are less lu- 
minous in X-rays than the Pleiades members, and suggest 
that this may be attributed to slightly older age, or lower 
rotation rates. 



1.1 Evolution of stellar angular momentum 

For a discussion of the context of t his work, the reader is re- 
ferred to our M34 publication in llrwin et"al] (|2006l ). The 
present survey probes a different open cluster of slightly 
younger age ('^ 150 Myr for NGC 2516, compared to 
~ 200 Myr for M34). Due to the small distance modulus 
and our use of a larger telescope, the NGC 2516 survey 
probes lower masses, covering a range of 0.1 < M/Mq < 0.7, 
compared to 0.25 < M/M© < 1.0 in M34, allowing the very 
low-mass stars to be examined in much greater detail. This 
will allow models of angular momentum evolution to be 
tested over a wider range of stellar mass, in particular to 
dete rmine whether models incorporati ng solid body rotation 
(e.g. lBouvier. Forestini fc Allainl|l997^ are able to reproduce 
the observed low-mass angular momentum evolution. At the 
lowest masses, core-envelope decoupling should not be rele- 
vant, since all the stars will be fuUy-convective and thus ro- 



tate approximately as solid bodies, allowing us to constrain 
directly the angular momentum loss law without needing 
to consider the details of internal angular momentum trans- 
port. At the age of NGC 25 16, stars with masses M < 0.4 M© 
should be fuUy-convective l|Chabrier fc Baraffelll997l ). 



1.2 The survey 

We have undertaken a photometric survey in NGC 2516 us- 
ing the CTIO 4m Blanco telescope and Mosaic-II imager. 
Our goals are two-fold: first, to study rotation periods for 
a sample of low-mass members, covering M spectral types, 
down to ~ 0.1 Mq, and second, to look for eclipsing bi- 
nary systems containing low-mass stars, to obtain dynam- 
ical mass measurements in conjunction with radial veloci- 
ties from follow-up spectroscopy. Such systems provide the 
most accurate determinations of fundamental stellar param- 
eters (in particular, masses) for input to models of stellar 
evolution, which are poorly constrained in this age range. 
We defer discussion of our eclipsing binary candidates to a 
later paper once we have obtained suitable follow-up spec- 
troscopy. 

These observations are part of a larger photometric 
monitoring survey of young open clusters over a range of 
ages and metalicities (the Mon itor project; iHodgkin et al.l 
l2006l and [Aigrain et ahllin press! ). 

The remainder of the paper is structured as follows; the 
observations and data reduction are described in !j2l and the 
colour magnitude diagram (CMD) of the cluster and candi- 
date membership selection are presented in ^ The method 
we use for obtaining ph otometric periods is summarised in 
El (see H rwin et al. I I2OO6I for a more detailed discussion) , and 
our results are summarised in Sj5] We discuss the implications 
of these results in SJSI and our conclusions are summarised 
in El 

In several of the following sections, we use mass and 
radius estimates for the cluster members. We note here that 
these were derived from the 150 Myr NextGen models, us- 
ing the /-band magnitudes, rather than V — I colour or V 
magnitude, for reasons discussed in i]3.1l 



2 OBSERVATIONS AND DATA REDUCTION 

Photometric monitoring observations were obtained using 
the 4m CTIO Blanco telescope, with the Mosaic-II im- 
ager, during two four-night observing runs, separated by 
~ 1 week, in 2006 February-March. This instrument pro- 
vides a field of view of ~ 36' x 36' (0.37 sq.deg), using a mo- 
saic of eight 2k x 4k pixel CCDs, at a scale of 0.27" /pix. 
NGC 2516 was observed for ~ 8 hours per night. 

In order to maximise the number of cluster members 
covered by our survey, and given the large angular extent of 
NGC 2516, we used three fields, with a single 75 s i-band 
exposure in each, observed by cycling around all three fields, 
to give a cadence of ~ 9 minutes, covering ~ 1 sq.deg of the 
cluster (illustrated in Figure [S]). Our fields were chosen to 
be outside the central few arcmin of the cluster to avoid 
contamination of the images from a number of very bright 
(y 5 — 8) stars close to the cluster centre. 

The observing conditions were superb, with exceptional 
seeing (< 0.8" on a few nights) and ~ 6 photometric nights. 



© 0000 RAS, MNRAS 000, 000-000 



Rotation of low-mass stars in NGC 2516 3 




Figure 2. DSS image of NGC 2516 covering ~ 1.8° X 1.5°, cen- 
tred on 07 57 25 -60 55 00 (J2000), showing the coverage of the 
pr esent survey (number ed 8-chip mosaic tiles), the XMM survey 
of iPillitteri et ah and the optical survey of I Jeffries et al.l 

||2001| '). 
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Figure 3. Plot of RMS scatter as a function of magnitude for 
the i-band observations of a single field in NGC 2516, for all un- 
blended objects with stellar morphological classifications. The di- 
agonal dashed line shows the expected RMS from Poisson noise in 
the object, the diagonal dot-dashed line shows the RMS from sky 
noise in the photometric aperture, and the dotted line shows an 
additional 1.5 mmag contribution added in quadrature to account 
for systematic effects. The solid line shows the overall predicted 
RMS, combining these contributions. 



We also obtained a number of longer V^-band exposures in 
photometric conditions (3x 300s in each field) to generate 
a colour magnitude diagram (CMD). Our observations are 
sufficient to give 1% or better pliotometric precision pre data 
point from saturation at i ~ 15 down to i ~ 19 (see Figure 
[3|, covering early to mid M spectral types at the age and 
distance of NGC 2516. 

For a full descr iption of ou r data red uction steps, the 
reader is referred to llrwin et al.l (|in pressi) . Briefly, we used 
the p ipeline for the INT wide-field survey (jlrwin fc Lewis! 
I2OOI ') for 2-D instrumental signature removal (crosstalk cor- 
rection, bias correction, flatfielding) and astrometric and 



photometric calibration. We then generated the master cat- 
alogue for each filter by stacking 20 of the frames taken 
in the best conditions (seeing, sky brightness and trans- 
parency) and running the source detection software on the 
stacked image. The resulting source positions were used to 
perform aperture photometry on all of the time-series im- 
ages. We achieved a per data point photometric precision of 
~ 2 — 4 mmag for the brightest objects, with RMS scatter 
< 1% for i< 19 (see Figure O. 

Our source detection software fiags as likely blends any 
objects detected as having overlapping isophotes. This in- 
formation is used, in conjunction with a morphological im- 
age classification fiag also generated by the pipeline software 
(|lrwin fc Lewisir2001 ) to allow us to identify non-stellar or 
blended objects in the time-series photometry. 

Photometric calibratio n of our data w as carried out us- 
ing regular observations of iLandoltl (|l992l ) equatorial stan- 
dard star fields in the usual way. 

Lightcurves were extracted from the data for ~ 100 000 
objects, 63 000 of which had stellar morphological classifica- 
tions, using o ur standard apertur e photometry techniques, 
described in ilrwin et al.l (|in pressi) . We fit a 2-D quadratic 
polynomial to the residuals in each frame (measured for each 
object as the difference between its magnitude on the frame 
in question and the median calculated across all frames) as 
a function of position, for each of the 8 CCDs separately. 
Subsequent removal of this function accounted for effects 
such as varying differential atmospheric extinction across 
each frame. Over a single CCD, the spatially-varying part 
of the correction remains small, typically ~ 0.02 mag peak- 
to-peak. The rea sons for using this tec hnique are discussed 
in more detail in llrwin et al.l l|in pressi) . 

For the production of deep CMDs, we stacked 20 i-band 
observations, taken in good seeing and photometric condi- 
tions. The limiting magnitudes, measured as the approxi- 
mate magnitude at which our catalogues are 50% complete 
(see i|3.2|l on these images were V ~ 24.3 and i ~ 22.7. 



3 SELECTION OF CANDIDATE LOW-MASS 
MEMBERS 

Catalogues of phot ometrically-selected candidate members 
were available from ljeffries et al] l|200ll ) for so lar mass down 
to ~ 0.2 M© and for 0.2 - 0.05 M© from iMoraux etHI 
(2005), but we elected to perform a new photometric selec- 
tion using V versus V ~ I CMDs from our data, to match 
exactly the field of view and magnitude range of our time- 
series photometry. Unfortunately, since we have avoided the 
very centre of the cluster in our field selection to reduce the 
effect of saturated bright stars on our CCD frames, the over- 
lap in sky coverage be tween the present s urv ey and many of 
the previous ones (e.g. lJeffries et al.ll200ll and lPillitteri et ahl 
2006) is relatively poor. 

3.1 The V versus V -I CMD 

Our CMD of NGC 2516 is shown in Figure H The V and 
i measurements were converted to the standard Johnson- 
Cousins photometric system using colour equations derived 
from our standard star observations: 



iv-i) 



{Vccd - iced) I 0.899 



(1) 
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Figure 4. V versus V - I CMD of NGC 2516 from stacked im- 
ages, for all objects with stellar morphological classification. The 
cluster main sequence is clearly visible on the right-hand side 
of the diagram. The boundaries of the region to select photo- 
metric candidate members are shown by the dashed lines (all 
objects between the dashed lines were selected). The reddening 
vector for Ay = 1.0 is shown at the right-hand side of the di- 
agram. The mass s cale is from the 150 Myr NextGen models 
llBaraffe et al.lll998h for M > 0.1 Mq, and the 120 Myr DUSTY 
models l lChabrier et al.ll2000l) for M < 0.1 Mq, using our em- 
pirical isochrone to convert the V magnitudes to / magnitudes, 
and subsequently obtaining the masses from these, due to known 
problems with the V magnitudes from the models (see i|3.1|l . The 
error bars at the left-hand side of the plot indicate the typical 
photometric error for an object on the cluster sequence. 



V = Kcd + 0.005 (V - /) (2) 
/ = iced - 0.096 (V - /) (3) 

Candidate cluster members were selected by defining an 
empirical main sequence 'by eye' to follow the clearly-visible 
cluster single-star sequence. The cuts were defined by mov- 
ing this line along a vector perpendicular to the cluster se- 
quence, by amounts k—aiV —I) and k~\-a{y —I) as measured 
along this vector, where aiV — /) is the photometric error 
in the V — I colour. The values of k used were —0.15 mag 
for the lower line and 0.5 mag for the upper line on the 
diagram, making the brighter region wider to avoid reject- 
ing binary and multiple systems, which are overluminous for 



16 18 20 22 24 26 

V magnitude 

Figure 5. Completeness in our source detection, measured as 
the fraction of simulated objects which were detected in each 
magnitude bin. The diagram is plotted in the Johnson V system 
assuming objects to lie on the empirically-derived cluster main 
sequence track of i|3.1l and Figure |4] The diagram has been ter- 
minated at V = 26 since we have not continued the selection of 
cluster members past this magnitude due to termination of the 
apparent main sequence in the CMD of Figure l4l 



their colour compared to single stars. 1685 candidate pho- 
tometric members were selected, over the full V magnitude 
range from V = 15.5 to 26, but the well-defined main se- 
quence appears to terminate at the hydrogen burning mass 
limit of M ~ 0.072 Mq, or V ~ 24, with a few candidate 
brown dwarfs found below this limit, but with high field 
contamination. 

We also consi dere d using the model i sochrones of 
iBaraffe et al] (|l998l ) and lChabrier et all (|2000l ') for selecting 
candidate members. The NextGen model isochrones were 
found to be unsuitable due to the known discrepancy be- 
tween these models and observations in the V — I colour 
for Teff <3700 K (corresponding here to V - I>2) This 
was examined in more detail by iBaraffe et al.l l|l998l ). and 
is due to a missing source of opacity at these temperatures, 
leading to overestimation of the 1/-band flux. Consequently, 
when we have used the NextGen isochrones to determine 
model masses and radii for our objects, the 7-band absolute 
magnitudes were used to perform the relevant look-up, since 
these are less susceptible to the missing source of opacity, 
and hence give more robust estimates. 



3.2 Completeness 

The completeness of our source detection was estimated by 
inserting simulated stars at random x,y positions into our 
images, drawing the stellar magnitudes from a uniform dis- 
tribution. Figure [S] shows the resulting plot of completeness 
as a function of V-band magnitude. The completeness for 
objects on the cluster sequence is > 90% down to the BD 
limit {V ~ 24). We suspect that the gradual drop in com- 
pleteness toward fainter magnitudes in Figure [5] results from 
a variable background in our CCD images caused by halos 
and scattered fight from the many bright stars in the cluster 
field, impeding detection of faint sources in certain regions 
of the frame. 
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Figure 6. Contamination, measured as the ratio of the calculated 
number of objects in each magnitude bin from the Galactic mod- 
els, to the number of objects detected and classified as candidate 
cluster members in that magnitude bin. 



3.3 Contamination 

In order to estimate the level of contamination in our cata- 
loRue , we used the Besangon Galactic models (|Robin et alj 
I2OO3I ) to generate a simulated catalogue of objects passing 
our selection criteria at the Galactic coordinates of NGC 
2516 (/ = 273.8°, h = -15.9°), covering the total FoV of 
~ 1.1 sq.deg (including gaps between detectors). We selected 
all objects over the apparent magnitude range 14 < 1/ < 26, 
giving 105 000 stars. The same selection process as above 
for the cluster members was then applied to find the con- 
taminant objects. A total of 651 simulated objects passed 
these membership selection criteria, giving an overall con- 
tamination level of ~ 39%. Figure |6] shows the contamina- 
tion as a function of V magnitude. We note that this figure 
is somewhat uncertain due to the need to use Galactic mod- 
els. Follow-up data (spectroscopy) will be required to make 
a more accurate estimate. We note in particular that the 
apparently very high contamination in the V = IQ — 17 bin, 
and possibly also the V = 2^ — 2Q bin, is probably the re- 
sult of the model poorly-fitting the observations, causing the 
contamination level to be overestimated. 



3.4 Binary fraction 

Since the cluster binary sequence is so clearly visible in 
the CMD of Figure U) we have computed a simple pho- 
tometric lower limit for the binary fraction in the cluster 
as a function of mass, using the empirical isochrone to se- 
lect stars with luminosities lying within a band of height 

2.5 logjQ 2 ~ 0.75 mag centred on the single star cluster se- 
quence, and the binary sequence (obtained by shifting the 
cluster sequence by this amount to greater luminosities), il- 
lustrated in Figure [7] This technique suffers from several 
sources of systematic error, including: (a) since the single 
star sequence is closer to the maximum of the distribution 
of field stars in the CMD, the single star counts are increased 
to a greater degree than the binary counts by field contam- 
ination, (b) binaries with mass ratios significantly different 
from 1 lie closer to the single star sequence, and therefore 
tend not to be selected as binaries, and (c) at the faint end of 
the CMD the binaries are brighter, and hence more readily 
detected than single stars of the same colour. 
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Figure 7. V versus V—I CMD, as Figure|4] showing the selection 
of photometric candidate single stars (grey) and binaries (black). 
The solid line shows the original empirical cluster main sequence, 
and the dashed lines show the boundaries used for selection, with 
y-band magnitudes V -f 1.25 logj^g 2, V — 1.25 logj^Q 2, and V — 
3.75 logj^Q 2, where V is the empirical main sequence magnitude. 
Typical photometric errors are indicated as for Figure |4] by the 
error bars at the left-hand side of the diagram. 



The binary fraction is shown as a function of mass in 
Figure [S] A total of 357 candidate binaries and 1065 can- 
didate single stars were selected, giving an overall binary 
fraction of 33.5 ± 2.1 % (using Poisson errors). 

We have been un able to use the more traditional I — K 
CMD approach (e.g. iPinfield et al] |2003| ) to determining 
limits on the binary fraction due to the lack of suitable near- 
IR data. 2MASS provides _ft'-band photometry, but the re- 
sulting I — K CMD did not show a clearly-visible binary 
sequence due to the large scatter caused by the photometric 
errors at the faint end of 2MASS. 

At this point, it is prudent to examine the sensitivity 
of our selection to the mass ratio q of the binaries. To do 
this, we have simulated binaries over a range of q from 
to 1.0, using the PMS models to compute the total J-band 
luminosity of the binary, and our empirical isochrone to con- 
vert the individual J-band magnitudes of the primary and 
secondary to V-band, and hence place the objects on the 
V — I CMD. Our selection was then applied to determine 
the minimum mass ratio q for the detected binaries, defined 
as 5 = M2/M1 (where Mi and M2 are the primary and 
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Figure 8. Photometric binary fraction as a fun ction of stel- 
lar m ass (computed using the NextGen models of iBaraffe et al.l 
Il998l) . 
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Figure 9. Minimum mass ratio q for detection of binaries as 
a function of the primary mass M\ using our CMD approach. 
Computed using the NextGen (for M > 0.1 Mq) and DUSTY 
jChabrier et al.ll2000l) models (for M < 0.1 Mq) to convert be- 
tween mass and 7-band absolute magnitude, and our empirical 
main sequence isochrone to convert between /-band and V-band 
magnitudes. The curve is not smooth owing to the coarse binning 
used in our empirical isochrone. 



secondary masses, respectively). In general we find that our 
selection is sensitive to smaller mass ratios at lower masses, 
which could explain the form of the binary fraction as a 
function of mass seen in Figure ID The figure indicates that 
we are sensitive to g>0.7 over the entire mass range, and 
at best, q > 0.3 at M ^ 0.2 M q. Comparing these values to 
those of iPinfield eFaLl (|2003l ), we find that the V — I selec- 
tion is slightly less sensitive to unequal mass ratio systems 
(e.g. in their 0.6— 1.0 M© bin the minimum detectable mass 
ratio was q — 0.5 whereas here it is g = 0.6). 

We note that since the completeness in Figure [S] is close 
to constant over the entire mass range covered by Figure |8l 
the apparent increase in binary fraction toward lower masses 
is probably not an incompleteness effect resulting from lower 
completeness for the single star sequence, compared to the 
more luminous binary sequence. However, our selection is 
affected by field contamination, with a background of field 
objects falling inside the selection bins on the CMD. This 
will tend to decrease the computed binary fraction, since 



there are more field objects at the single star side of the 
selection region (see Figure |4]). It is difficult to correct for 
this effect until we have obtained follow-up spectroscopy. 

Our res ults in V — I are in good general agreement 
with those of lPinfield et al.1 (120031 ). although it is difficult to 
distinguish whether our data are consistent with any change 
in the binary fraction as a function of mass, given the small 
number statistics (and hence large Poisson errors in Figure 
[5]) and better sensitivity to small q with decreasing mass 
(leading presumably to more detections). 

Photometric selection of binaries remains, to our knowl- 
edge, relatively untested, without detailed evaluation of the 
sensitivity. Therefore, a more detailed and reliable study of 
binarity in NGC 2516 requires follow-up spectroscopy. 



4 PERIOD DETECTION 
4.1 Method 

The method we use f or detection of pe riodic variables is 
described in detail in llrwin et al.l l|2006l ), and we provide 
only a brief summary here. The method uses least-squares 
fitting of sine curves to the time series m{t) (in magnitudes) 
for all candidate cluster members, using the form: 



m{t) — TUdc + a sin{ujt + t 



(4) 



where rUdc, Q (the amplitude) and <j!> (the phase) are free 
parameters at each value of uj over an equally-spaced grid of 
frequencies, corresponding to periods from 0.005 — 50 days 
for the present data-set. 

Periodic variable lightcurves were selected by evaluating 
the change in reduced x^'- 



Axt 



smooth 



> 0.4 



(5) 



where Xi/ is the reduced x of the original lightcurve with 
respect to a constant model, and xt smooth is the reduced 
X^ of the lightcurve with the smoothed, phase-folded ver- 
sion subtracted. This threshold was used for the M34 data 
and appears to work well here too, carefully checked by ex- 
amining all the lightcurves for two of the detectors, chosen 
randomly. A total of 1011 objects were selected by this au- 
tomated part of the procedure. 

The selected lightcurves were examined by eye, to de- 
fine the final sample of periodic variables. A total of 362 
lightcurves were selected, with the remainder appearing non- 
variable or too ambiguous to be included. 



4.2 Simulations 

Monte Carlo simulati ons were performed following the 
method detailed in Ir win et al.l ([2006), injecting simulated 
signals of 2% amplitude and periods chosen following a uni- 
form distribution on logjQ period from 0.1 to 20 days, into 
lightcurves covering a uniform distribution in mass, from 0.7 
to 0.1 Mq. a total of 2271 objects were simulated. 

The results of the simulations are shown in Figure [TU] as 
greyscale diagrams of completeness, reliability and contam- 
ination as a function of period and stellar mass. Broadly, 
our period detections are close to 100% complete down 
to 0.2 M0, with remarkably little period dependence. The 
slight incompleteness (approximately 63% completeness for 
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Figure 11. Detected period as a function of actual (input) period 
for our simulations. Objects plotted with crosses had fractional 
period error < 10%, open circles > 10%. The straight lines rep- 
resent equal input and output periods, and factors of 2, 3, 1/2 
and 1/3. The curved lines are the loci of the ±1 day~^ aliases 
resulting from gaps during the day. The majority of the points 
fall on (or close to) the line of equal periods. 

0.6 < M/Mq < 0.7) seen in the highest mass bin is a result 
of saturation. Figure [TT] shows a comparison of the detected 
periods with real periods for our simulated objects, indicat- 
ing remarkably high reliability, especia lly compared to the 
earlier M34 data-set (|lrwin et al.ll2006l ). resulting from bet- 
ter sampling. 



4.3 Detection rate and reliability 

The locations of our detected periodic variable candidate 
cluster members on a, V versus V - I CMD of NGC 2516 
are shown in Figure \n\ The diagram indicates that the ma- 
jority of the detections lie on the single-star cluster main 
sequence, as would be expected for rotation in cluster stars 
as opposed to, say, eclipsing binaries. The measured pho- 
tometric binary fraction (see i|3.4|l for the periodic variable 
stars is 26.3 ± 3.5% (using Poisson errors). This value is 
slightly lower than the value for all cluster members from 
i|3.4l but only at the 1.8a level. Given the systematic un- 
certainties in the binary fractions (which are not accounted 
for by the Poisson counting errors), this is probably not a 
significant result. 

Figure [13] shows the fraction of cluster members with 
detected periods as a function of V magnitude. The decay- 
ing parts of the histogram at the bright and faint ends may 
be caused by the increased field contamination here (see 
Figure O, since we expect field objects on average show less 
rotational modulation than cluster objects, and/or incom- 
pleteness effects resulting from saturation for V < 17, and 
for V > 21, the gradual increase in the minimum amplitude 
of variations we can detect (corresponding to the reduction 
in sensitivity moving to fainter stars, see Figure [Sjl. 

In the M34 survey of llrwin et all (|200& ) , we commented 
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> 




Figure 12. Magnified V versus V - I CMD of NGC 2516, 
for objects with stellar morphological classification, as Figure |4l 
showing all 362 candidate cluster members with detected peri- 
ods (black points). The dashed lines show the cuts used to select 
candidate cluster members (sec i|3.1|l . 
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Figure 13. Fraction of candidate cluster members detected as 
periodic variables, plotted as a function of magnitude. This dis- 
tribution has not been corrected for incompleteness in the period 
detections, which are close to 100% complete for 17 < V < 21 
(shown by the vertical dashed lines). 
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Figure 10. Results of the simulations for 0.02 mag amplitude, plotted as a function of mass (left) and period (right). The simulated 
region covered 0.1 < M/M© < 0.7 in order to be consistent with the NGC 2516 sample. Top panels: completeness as a function of real 
(input) period. Centre panels: Reliability of period determination, plotted as the fraction of objects with a given true period, detected 
with the correct period (defined as differing by < 20% from the true period). Bottom panels: Contamination, plotted as the fraction 
of objects with a given detected period, having a true period differing by > 20% from the detected value. 



on a possible increase in the fraction of photometric variables 
from K to M spectral types, subject to a large uncertainty 
due to small number statistics. Since the mass range covered 
by the present sample is different, direct comparison is diffi- 
cult, and there is no clear evidence for such a trend in NGC 
2516. We further note that the decline observed in M34 at 
faint magnit udes is probably i ndeed due to incompleteness 
as stated by llrwin" et al.l (|2006l ). since a similar trend is also 
seen in the NGC 2516 data-set, displaced accordingly for the 
differences in sensitivity and distance modulus. 

We note that lower field contamination (as evalu- 
ated in H3.3P is expected in the rotation sample than in 
the full candidate membership sample. Typi cal field pop- 
ulation ages for the young disc of 3 Gyr l|Mayoij Il974l . 
iMeusinger. Stecklum fc Reimannlll991^ imply slower rota- 
tion rates by factors of a few than cluster members, re- 
duced activity, implying smaller asymmetric components of 
the spot coverage and hence lower photometric amplitudes 
(which probably render many of them undetectable). 

The properties of all our rotation candidates are listed 
in Table [T] 

4.4 Non-periodic objects 

The population of objects rejected by the period detec- 
tion procedure described in S|4] was examined, finding that 
the most variable population of these lightcurves (which 



might correspond to non-periodic or semi-periodic variabil- 
ity) was contaminated by a small number of lightcurves 
(~ 50) exhibiting various uncorrected systematic effects, 
mostly seeing-correlated variations due to image blending, ft 
is therefore difficult to quantify the level of non-periodic or 
semi-periodic variability in NGC 2516 from our data. Qual- 
itatively however, there appear to be very few of these vari- 
ables, and examining the lightcurves indicated only ~ 10 
obvious cases, some of which resembled eclipses (either plan- 
etary transits or eclipsing binaries), and will be the subject 
of a later Monitor project paper. 

5 RESULTS 

5.1 NGC 2516 rotation periods 

Plots of period as a function oiV — I colour and mass for the 
objects photometrically selected as possible cluster mem- 
bers are shown in Figure 1141 These diagrams show a strik- 
ing correlation between stellar mass (or spectral type) and 
the longest rotation period seen at that mass, with a clear 
lack of slow rotators at very low masses. This trend is also 
followed by the majority of the rotators, with only a tail of 
faster rotators to ~ 0.25 day periods. Furthermore, very few 
objects were found rotating faster than this, implying a hard 
lower limit to the observed rotation periods at 0.25 days. 
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21.25 
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0.649 


0.013 


0.20 


0.24 


N2516-1-1-881 


07 56 47.57 
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0.010 


0.55 


0.51 


N2516-1-1-958 


07 56 49.99 


-61 32 19.9 


18.48 


15.99 


6.291 


0.015 


0.49 


0.45 


N2516-1-1-1470 


07 57 08.92 


-61 29 18.6 


17.67 


15.49 


8.803 


0.016 


0.56 


0.52 



Table 1. Properties of our 362 rotation candidates. The period P in days, t-band a mplitude ai (units of magnitudes, in the instrumental 
bandpass), interpolated mass and radius (from the models of ^B araffe et a l.^ 1998, derived using the / magnitudes) are given (where 
available). Our identifiers are formed using a simple scheme of the cluster name, field number, CCD number and a running count of stars 
in each CCD, concatenated with dashes. The full table is available in the electronic edition. 
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Figure 14. Plots of rotation period as a function of dereddened 
V — I colour (top), and mass (bottom) for NGC 2516, deriving the 
masses using the 15 Myr NextGen mass-magnitude relations of 
iBaraffe et"aj1 l ll99^ and our measured 7-band magnitudes. In the 
lower diagram, the greyscales show the completeness for 0.02 mag 
periodic v ariations from our simulations, and X-ray sources from 
IPiUitteri et al. (20061 ) are overlaid with open circles. 



Could the apparent morphology in Figure [14] be ex- 
plained by sample biases? The simulations of M.2\ suggest 
that this is unlikely, since we are sensitive to shorter periods 
than the 0.2 — 0.3 day 'limit', and the slight bias toward 
detection of shorter periods at low mass is not sufficient to 
explain the observations for the slow rotators. Furthermore, 
Figure [15] indicates that the lack of sensitivity to low am- 
plitudes at low masses does not appear to introduce any 
systematic changes in the detected periods. 

Figures[l4]and[T5]show an apparent lack of objects with 
masses 0.2 < M/Mq < 0.4 and rotation periods > 2 days, 



--0.1 1 10 

Period (doys) 

Figure 15. Plot of amplitude as a function of period for NGC 
2516 in two mass bins: 0.4 < M/Mq < 0.7 (top) and 0.2 < 
M/Mq < 0.4 (bottom). 

a region of the diagram where our survey is sensitive down 
to amplitudes of 0.01 mag (see also H4.2p . There is also an 
apparent lack of large-amplitude objects in this region of the 
diagram, although this is not statistically significant (apply- 
ing a Kolmogorov-Smirnov test between the two mass bins 
indicated a probability of 0.25 that the distributions were 
drawn from the same parent population). 

5.1.1 Period distributions 

In order to quantify the morphology of Figure 1141 we have 
used histograms of the rotation period distributions in two 
broad mass bins, 0.4 < M/Mq < 0.7 and M < 0.4 Mq, 
shown in Figure [161 We have attempted to correct the distri- 
butions for the efi'ects of incompleteness and (un)reliability 
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using th e simulations descri bed in t]4.2l following the method 
used in llrwin et all (|2006l ). The results of doing this are 
shown in the solid histograms in Figure 1161 and the raw 
period distributions in the dashed histograms. 

The period distributions in the two mass bins of Fig- 
ure [16] show clear differences, with the low-mass stars {M < 
0.4 M0) showing a strongly peaked rotational period distri- 
bution, with a maximum at ~ 0.6 — 0.7 days, whereas the 
higher-mass stars (0.4 < A'I/Mq < 0.7) show a broader dis- 
tribution. We applied a two-sided Kolmogorov-Smirnov test 
to the corrected distributions to confirm the statistical sig- 
nificance of this result, finding a probability of 4 x lO"^'' that 
the distributions were drawn from the same parent popula- 
tion. 

The implication of this result is that the observed mor- 
phology in Figure 1141 and in particular the increase of the 
longest observed rotation period as a function of increasing 
mass, a trend followed also by the bulk of the rotators, is 
real and statistically significant. 



5.1.2 Low-mass slow rotators 

A group of ^ 4 outliers are visible e.g. in the upper 
panel of Figure 1141 as unusually slow rotators at low mass 
{(V ~ I)^ 2.5). These objects also give rise to the apparent 
'trend' seen in Figure [15] of smaller amplitudes for periods 
P>2 days, for the lower mass bin (0.2 < M/Mq < 0.4). 

Careful examination of the lightcurves indicates that 
these objects are all convincing detections, and lie on, 
or close to, the cluster sequence. Obtaining optical spec- 
troscopy for these objects would appear to be warranted, 
to constrain their membership of the cluster, and determine 
spectral types, to evaluate their reddening. 



5.1.3 Rapid rotators 

We have examined the periods of our fastest-rotating stars, 
to check if they are rotating close to their break-up velocity. 
The critical period Pcrit for break-up is given approximately 
by: 



0.116 days 



N3/2 



(M/M0)i/2 



(6) 



where R and M are the st ellar radius and mass respec- 
tively (e.g . 'Herbst et al]l2002l ). Using the NextGen models of 
iBaraffe et al. (1993)7 the object rotating closest to breakup 
is N2516- 2-8-624, with Pcrit/P = 0.42. However, the major- 
ity of our objects are rotating at much lower fractions of 
their break-up velocity. 

Figure [TS] indicates that all but three of our fastest 
rotators have P > 0.2 days. The outliers are N2516-1-4- 
3373 (P = 0.148 days), N2516- 2-8-566 (P = 0.146 days) 
and N2516- 2-8-624 (P = 0.087 days). Two of these objects, 
N2516- 1-4-3373 and N2516- 2-8-624, clearly lie on the binary 
sequence in Figure|4] and N2516- 2-8-566 may also lie on the 
binary sequence, so we suspect that these stars might be 
tidally locked binaries. They are therefore good candidates 
for spectroscopic follow-up. 

We note that there are three outliers visible in Figure 
[14] with {V — 1)0 < 1.3, and periods close to 1 day. Sim ilar 
objects were seen in our M34 study (,Irwin et al.n2()06l ). At 



these masses, field contamination in our survey is signifi- 
cant (~ 60%), so we suspect that they may be contaminat- 
ing field objects, and possibly binaries (in M34, several of 
the corresponding objects were also detected in X-rays, it is 
difficult to say if this holds in NGC 2516 due to the incom- 
pl ete overlap in spatia l coverage with the X-ray catalogue 
of Pi llitteri et al1l2006l ). Our survey is incomplete at these 
(high) masses due to saturation of the detectors, which var- 
ied between CCDs in the mosaic, giving rise to the lack of 
other objects detected in this range of {V — I)o. 



5.2 Comparison with other data-sets 

5.2.1 Period versus mass 



Figure [T7] shows a diagram of rotation perio d as a func- 
tion of stellar mass for the O NC (1 ± 1 Myr; iHillenbrandl 



1 199711, NGC 2 264 f2 - 4 Mvr: lPark et al.ll2000l). NGC 236 2 



(~ 5 ± 1 Myr; iMoitinho eraLll200ll, iBalona fc La"nevlll996l '). 
the Pleiades (~ 100 Myr: iMevnet et al.lll993l ). NGC 2516 
and M34 (~ 200 Mvr: lJones fc Prosserlll996l '). Data sources 
for each cluster are indicated in the figure caption. From this 
Figure, it appears that the objects we observe in NGC 2516 
are bounded by lines of P ~ constant (or a low power of Al) 
for the fastest rotators, and P oc for the slow rotators 
in the mass range 0.2 < M/Mq < 0.7. 

Furthermore, comparing the NGC 2516 and M34 results 
in Figure [T7] indicates that the slowest rotators for M > 
0.6 M0 have P ~ constant as a function of mass, whereas 
for A'l < 0.6 Mq, the period increases as a function of mass. 
A similar trend is visible in the Pleiades v sin i data (noting 
that a number of the slowest rotators here only have v sin i 
lower limits, denoted by double-headed arrows). 

The diagram clearly shows a gradual evolutionary se- 
quence, from a relatively fiat mass-dependence of the rota- 
tion periods in the ONC (~ 1 Myr), to a sloping relation 
in NGC 2362 (~ 5 Myr), and the emergence of the break 
between a flat distribution for M > 0.6 Mq and strongly 
sloping distribution at lower masses, for the Pleiades, NGC 
2516 (~ 150 Myr) and M34 (~ 200 Myr). We suggest this 
apparent change in slope at M ~ 0.6 M© may correspond to 
a transition to convectively-dominated stellar interiors mov- 
ing to lower masses (but see ij6.2|l . 

5.2.2 V sin i measurements 



iTerndrup et al.l (|2002l ) presented v sin i measurements for G 
and early-K dwarfs in NGC 2516. Dividing their sample into 
branches of fast and slow rotation for {B — V)o > 0.6 and 
counting objects with vsini < 15 km s~^ as slow rotators, 
gave 7 fast rotators, with typical usini ~ 20 — 50 km s~^, 
and 17 slow rotators with ?;sini< 10 km s~^, i.e. 30% are 
fast rotators. The threshold of vsini — 15 km s~^ corre- 
sponds to P = 2.2 sini days at 0.7 M0, or (P) = 1.72 days, 
using (sini) — n/A (averaged over all inclinations i). Apply- 
ing this threshold to the 0.5 — 0.7 Mq bin of Figure [Til we 
counted 69 slow rotators, and 40 fast rotators, i.e. 37% fast 
rotators (and clearly from Figure [14] this fraction increases 
considerably for lower masses, to nearly 100% at 0.2 Mq). 
Furthermore, for {B — V)o < 0.6 there are no clearl y-visible 
fast rotators in Figure 5 of iTerndrup et all ||2002D . These 
results clearly indicate that the fraction of objects rotating 



© 0000 RAS, MNRAS 000, 000-000 



Rotation of low-mass stars in NGC 2516 11 



0.4 < M/Mp, <= 0,7 



0.4 < M/Mq <= 0.7 



Period (days) 
M < 0.4 M„ 



20 



Period (days) 
M < 0.4 M„ 



±11 



Period (days) 



Period (days) 



Figure 16. Period distributions for objects classified as possible photometric members, in two mass bins: 0.4 < M/M© < 0.7 (upper 
row, corresponding rougiily to early-M spectral types) and M < 0.4 (lower row, late-M). The left-hand panels show the distributions 
plotted in linear period, and the right-hand panels show the same distributions plotted in logjg period. The dashed lines show the measured 
period distributions, and the solid lines show the results of attempting to correct for incompleteness and reliability, as described in the 
text. 



faster than the chosen threshold increases as a function of 
mass, lending support to the rotation period results we have 
presented. 



5.2.3 Comparison to the field: main sequence spin-down 



Comp aring our results with the field survey of lDelfosse et alj 
l| 19981 ) indicates that the spin-down timescale of M- 
dwarfs increases moving to later spectral types (or lower 
masses). The fastest rotators observed by t hese authors 
in the young disc population (age ~ 3 Gyr; iMavoj Il974l . 
iMeusinger. Stecklum fc Reimannll99ll ) were of ~ M4 spec- 
tral type, with P/ sini ^ 12 hours. Given that we expect 
the objects to spin-down on the main sequence, their origin 
must be in the most rapid rotators of the present sample 
(~ 0.25 days), suggestin g a factor of only ~ 2 spin-down. 
This is slower than the SkumanichI (| 19721 ) t^^^ prediction 
of a factor of ^ 4.5, suggestive of a very long spin-down 
timescale for these late-type dwarfs compared to solar-type 
dwarfs. 



6 DISCUSSION 

6.1 Link to angular momentum 

It is instructive to examine the implications of the relations 
determined in >i5.2.1l for the stellar angular momentum J. 
We define 



J = luj 

27rfc2 Mi?2 



(7) 

(8) 



where / = k^MR^ is the moment of inertia of a star of mass 
M and radius R, k is the radius of gyration (e.g. fc^ = 2/5 for 
a sphere of uniform density), and u = 2-k jP is the rotational 
angular velocity. 

For a fuUy-convective star (Af <0.4 M©), k is approx- 
imately independent of mass, and therefore the specific an- 
gular momentum j — J/M oc R? /P. Furthermore, on the 
ZAMS we can approximate M oc R at these masse^J Thus, 



Fitting th e ma ss-radius relation from the models of 



iBarafFe et al] l| l99^ over the range 0.1 < M/Mq < 0.7 with 
a simple linear model gave R/Rq = O.SAZ/Mq -|- 0.075, which is 
a better approximation than the simple proportionality. 
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Figure 17. Rotation period as a function of stellar mass for 
(top to bottom): ONC, NGC 2264, NGC 2362, the Pleiades 
(rotation periods and vsini data on separate plots), NGC 
2516 and M34. Lower and upper limits (from vsini data) 
are marked with arrows. The masse s were taken from the 
NextGen mass-magnitude relations ijBaraffe et al.l Il998l) at 
the a ppropriate ages. The ONC data are from ]Herbsret_ajJ 
ll2002l'). For NGC 2264 we used the data of iLamm et al.l l|2005l 'l 
and iMakidon et al.l | |2004|) . The NGC 2362 data are from the 
Monitor project, to be published in Hodgkin et al. (in prep). The 
Pleia des rotation period data are a compilatio n of the results 
from IVan Locuwcn. AlDhcnaar fe^e^d dlQSTlV Istauffcr ct al. 

IQSTllTlMagnitskii (1987 

1993bl ). |Prosser et al.ll|l995h . iKrishnamurthi et al.l|l998l'l f taken 



from the open c l uster d atabase), lTerndrup''et''alT ill 9991) and 
IScholz fc Eisloffell ll2q04bl'). The Ple i ades v sin i data are a compi 
lation of results fromlStauffer et "al] lll984l').IStauffer &: HartmannI 



1987]') . ISoderblom et al.l lll993l]. Ijones. Fischer fc Stauffei 
19961) . iQueloz et al.l lll998^ and iTerndrup et al.l (|2000|V The 



M34 data are taken from llrwin et alJ l|2006h . In the NGC 2516 




M / Mg 

Figure 18. Stellar angular momentum J as a function of mass 
for the NGC 2516 rotation period sample. The lines show J = 
constant, J on M and J oc . The error bars show the median 
of the distribution binned in 0.1 dex bins of log Af. 



recalling the simple description of the morphology of the pe- 
riod versus mass plots in terms of power law relations from 
>i5.2.1l and Figure 1171 P = constant implies j oc , and 
P cx implies j oc . These correspond to J oc 
and J = constant, respectively. A plot of J as a function of 
stellar mass for the NGC 2516 sample is shown in Figure [THl 
with these lines overlaid. 

We computed the median rotation period as a function 
of mass for the NGC 2516 data, finding it was consistent 
with Pmod oc M^, a relation inbetween those of the upper 
and lower envelopes in Figure [17] This implies j — constant 
as a fu nction of mass, a result consistent wit h the conclu- 
sions oflHer bst. Bailer- Jones fc Mundtl lj200ll l in the ONC. 
The lack of any change in the median relation for stars still 
undergoing contraction on the PMS indicates that any mass- 
dependent angular momentum losses are not important in 
determining the evolution of the median over this age range. 
However, this is not the case in detail or once the stars reach 
the ZAMS, in particular the shape of the distribution even 
in this mass range does evolve in a mass-dependent fashion. 

It is interesting to co mpare the present observ ations 
with the Pleiades survey of IScholz fc Eisloffcl' (2004b'). The 
mass-dependence of the angular momentum evolution we 
have derived for NGC 2516 appears to be slightly incon- 
sistent with their conclusions, where they found P oc M. 
Examining the Pleiades diagram in Figure [17] indicates that 
their result can most likely be attributed to the lack of data 
for 0.25 < M/Mq < 0.5. Otherwise the results of the sur- 
veys are in good agreement, as expected due to the similar 
ages of the two clusters. 

The earlier conclusions of IScholz fc Eisloffel l|2004ah in 
a Ori (age ~ 3 Myr) are interesting, since they also found 
P oc M. At these ages, the mass-radius relation is better-fit 
by R/Rq oc (M/Mq)'^''^, so this implies P oc i?^, or con- 
stant specific angular momentum, the same as our median 
relation. 
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6.2 Simple models 

6.2.1 Description 

We have developed a simple model scheme, to attempt to 
find an e xplanation for the ob servat ions, fo llowing the ap- 
proach of iBouvier et all (119971 ) and lAllainI ((l998). In this 
work, the relation between M, R and fc as a function of 
time for the radiative core (where present) and convective 
envelope was taken from the Lyon group stellar models, gen- 
erously supplied to us by I. Baraffe (private communication) 
for 0.1,0.2,0.3,0.4,0.6,0.8 and 1.0 Mq. We have interpo- 
lated the values where necessary for intermediate masses. 

The overall angular momentum loss rate was split into 
two components: losses due to stellar winds, assuming a loss 
law with saturation at a critical angular velocity tJsat (al- 
lowed to vary as a function of mass), and losses due to disc 
locking, which to a good approximation maintains a con- 
stant angular velocity until the circumstellar disc dissipates. 

We assume an angular m omentum loss l aw for stellar 
winds of the form suggested bv lKawalej (|l988t ). modified to 
include saturation of the ang ular momentum losses above 
a critical angular velocity Wsat . iBames. fc Sofial l| 19961 ) argue 
that saturation at high velocity is required to account for fast 
rotators on the ZAMS, and may be physically motivated by 
saturation of the magnetic dynamo activity giving rise to 
the angular momentum losses. This is supported by obser- 
vations of several diagnostics of magnetic activity, incl uding 
saturation of chromospheric emission (e.g. IVilhu|[l984l ). 

The adopted angular momentum loss law has the fol- 
lowing form; 



(-) 

\ dt J wind 



(4)) 

/ \ 1/2 / \ -1/2 



where u is the angular velocity of the star, R is its radius, 
and M its mass, if is a constant, chosen to give the correct 
value of the solar angular momentum for a 1 Mq star at 
the age of the Sun. The valu e K = 2.7 x 10*^ g cm^ s was 
determined bv lBouvier et all |l993) to satisfy this require- 
ment, with LJsat = 14 loq. This particular form of angular 
momentum loss law is found to reproduce the time depen- 
dence of rotation between the Hyades age (625 Myr) and 
the Sun (4.57 Gyr). 

The value of uisat is found to depend on stellar mass 
(e.g. iKrishnamurthi et~al]ll997l , iBarnes fc Sofialll996l ), and 
is typically assumed to be inversely proportional to the con- 
vective overturn timescale r: 

i^sa.t = t^sat,© — (10) 
T 

where the quantities subsc ripted are those for the sun 
l|Krishnamurthi et al.lll997|). The values of r were taken from 
the 200 Myr model of iKim fc Demarguel (Il996l ) for M > 
0.5 Mm. ISills. Pinsonneault fc Terndrud l|2000l ) found that 
simply linearly extrapolating these models to lower masses 
did not work well, and derived a set of empirical values to 
fit the distribution seen in the Hyades at these masses. We 
have also found this to be the case for t he NGC 2516 data, 
and adopt the values of lSills et al.](|200Gl ). which give a much 
better fit to the data. 

We consider two classes of models: solid body rotation, 
where the star is considered to rotate as a single solid body, 



implying that the core and envelope must have the same an- 
gular velocity, and differential rotation, where the core and 
envelope are allowed to decouple, with angular momentum 
losses from the envelope, and angular momentum transfer 
on a timescale Tc between the core and envelope. 

In order to account for disc locking, we introduce a disc 
locking time rdisc, such that for t < rdisc, the star is locked 
to the circumstellar disc, and maintains a constant angular 
velocity in the convective envelope, equal to the initial value. 
For t > Tdisc, the angular velocity evolves as follows: 

For solid body rotation, we can simply differentiate the 
expression J — lu where I, J and u} are the quantities 
for the entire star, with respect to time. The only angular 
momentum loss is due to stellar winds, so: 



(11) 



duj Lo / dj\ Lu dl 

dt J \ dt J wind I dt 

For differential rotation, the core and envelope must be 
considered separately. F or the envelope, two extra terms are 
added (see lAllair]|l998l for a full derivation), describing the 
coupling of angular momentum from the core (first term), 
and the loss of mass to the core as the star evolves (second 
term): 



dcjc. 



dt 



1 AJ 2 R: 

-^conv '^"c 3 /ci 



2 

rad 



dM,: 



conv 
^conv 

fdJ, 



dt 



dt ) wind 



dt 



where we assume: 

iconvt/rad irad<iconv 



(12) 



-^rad -^conv 

as suggested by iMacGregorl (|l99l[ ). 
(9) The corresponding equation for the evolution of the core 
angular velocity is: 



1 AJ 2 Ri 



dAL, 



<^rad din 



dt 



T r '^conv , 

J rad '^c J rad f^t' 



Jrad dt ^ 

where the first two terms again represent coupling of angular 
momentum to the envelope, and mass gained by the core 
(from the convective envelope), respectively. The final term 
accounts for the change in moment of inertia of the core as 
the star evolves. 



6.2.2 Evolution from NGC 2362 

We first attempt to evolve the observed rotation rates from 
our survey in NGC 2362 (Hodgkin et al., in prep) forward in 
time to reproduce the available rotation period data in NGC 
2516. In order to do this, we consider the slowest rotators, 
characterised by the 25th percentile (the lower quartile) of 
the distribution of observed angular velocities, uo, and the 
fastest rotators, characterised by the 90th percentile of the 
distribution. The discrepancy in measures here is because 
the 10 percentile was found to be somewhat unstable in the 
presence of small number statistics for slow rotators in some 
of the clusters under consideration. Nevertheless, the results 
for this quantity are compatible with the conclusions we 
draw for the 25 percentile for those clusters with sufficiently 
good statistics. 

We fit solid body and differentially rotating models to 
the data, allowing the parameters ojsat, Tdisc and Tc to vary 
to fit the slowest rotators (25 percentile). 
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Figure [19] shows our results, where we performed the 
analysis in a series of mass bins, in order to resolve the mass 
dependence of the evolution. These were: 0.9 < M/Mq < 
1.1, 0.7 < M/Mq < 0.9, 0.5 < M/Mq < 0.7, 0.35 < 
M/Mq < 0.5, and 0.2 < M/Mq < 0.35 (chosen empirically 
to encompass the changes in behaviour seen in the models, 
while retaining reasonable statistics). The models were cal- 
culated for single masses roughly at the centre of these bins, 
of 1.0, 0.8, 0.6, 0.42 and 0.28 M©. 

The solid body models shown in Figure [19] show good 
overall agreement with the upper 90 percentiles for the NGC 
2516 data, in all the mass bins. For the 0.9 < M/Mq < 1.1 
bin, the M34 upper 90 percentile is at much greater lj than 
the model, but we suspect that this result is not significant 
due to small number statistics in the M34 sample at these 
masses: the discrepancy i s due to th e pres ence of 3 data 
points, and as discussed in llrwin et al.l l|2006l) . these may be 
contaminants. For 0.7 < M/Mq < 0.9 the upper 90 per- 
centile in the Hyades appears to lie at much greater uj than 
the model, but again this result is strongly affected by small 
number statistics. In order to resolve these issues, improved 
statistics at these masses and ages will be required. We note 
that the differentially rotating models with the assumed val- 
ues of Tc give a poor fit to the fast rotators (as characterised 
by the 90 percentiles), especially on the ZAMS. 

For the lower 25 percentiles, the solid body models give 
a poor fit to the evolution from NGC 2362 to older ages, 
which becomes particularly apparent after ~ 50 Myr to the 
age of the Hyades in the top two mass bins. Introducing 
core-envelope decoupling with Tc = 1.0 Gyr provides a much 
better fit to the observations, and produces curves in good 
agreement with all the presently available data. We note 
however that we have not constrained these to fit the sun, 
which requires a modification of the constant K in Eq. ([9]). 

An interesting implication of the results in Figure [T9l is 
that we cannot find a model in this scheme that fits both 
the slow and fast rotators simultaneously: the values of Tc 
required in the models are different (solid body models cor- 
respond to Tc ~ 0). The result that the fast rotators seem to 
behave more like solid bodies may indicate that the redistri- 
bution of angular momentum at the core/envelope boundary 
is more efficient in fast rotators. 

We note that there are relatively few constraints on the 
evolution for masses < 0.7 Mq at the intermediate ages (5 
to 150 Myr), and at old ages (> 200 Myr). We have recently 
obtai ned data in NGC 2547 (~ 38.5 Myr; iNavlor fc Jeffries! 
I2OO6I ) at these masses, which should allow us to address the 
first of these problems in a forthcoming publication. 

6.2.3 Constraints on initial conditions 

We explore here an altern ative meth o d for analysis of the 
NGC 2516 data, foUowing ISiUs et all l|2000l ). By assuming 
that all stars sha re the sam e initial rotation rate, here 10 
days, as used by ISills et al.1 ()2000[ ). we can constrain the 
disc lifetimes required to reproduce the observed envelope in 
diagrams of o) as a function of stellar mass. In particular, the 
observed lower envelope as a function of mass (the slowest 
rotators seen) must result from the longest disc lifetimes, 
and hence we can derive the required maximum disc lifetime 
as a function of mass. Figure [20] shows the results of this 
analysis for NGC 2516. 




1 0.8 0.6 0.4 0.2 

M / M„ 




0.6 0.4 0.2 



M/Mq 

Figure 20. Top and centre: NGC 2516 dat a (points), plus 
lower limits derived from the v sin i values of iTerndrup et al.l 
||2002| ) (arrows), with models overlaid (top; solid body model, cen- 
tre: differentially rotating model) for tdisc = 0, 0.3, 1, 3, 10, 30 
and 100 Myr (top to bottom). The initial period Pq was 10 days. 
The straight line shows the approximate lower envelope used to 
derive the disc lifetime as a function of mass. Bottom: derived 
disc lifetime tdisc as a function of mass, for the solid body models 
(filled circles and solid line), and differentially rotating models 
(o pen circles and d ashed line). The models used the LJsat values 
of lsills et ai] ll2000t) (these were adjusted to provide a good fit to 
the Hyades data). 



The results in Figure [20l indicate that the disc lifetime 
must decrease strongly as a function of decreasing mass in 
order to reproduce the NGC 2516 data with a single ini- 
tial rotation period. Lifetimes of > 10 Myr are required for 
M > 0.3 Mq, a resul t at odds with other measurements 
(e.g. iHaisch et al.ll200ll : iBouwman et al ] |2006l and references 
therein), which indicate values of ~ 10 Myr at solar mass. 
Therefore it seems unlikely that the envelope shape we ob- 
serve is due purely to varying disc lifetimes as a function of 
mass. 

An alternative explanation (the other extreme) is that 
the initial period Pq varies as a function of mass, but the disc 
lifetime t^isc remains fixed. Figure [2l1 shows the results of do- 
ing this for the NGC 2516 data, assuming tdisc = 10 Myr, 
indicating again that a rather wide range of initial condi- 
tions appear to be required to explain the observed envelope 
purely by variations of Po- Furthermore, the very slow rota- 
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Figure 19. Rotational angular velocity lo plotted as a function of time in five mass bins: 0.9 < M/Mq < 1.1, 0.7 < M/Mq < 0.9, 
0.5 < A//Mq < 0.7, 0.35 < M/Mq < 0.5, and 0.2 < M/Mq < 0.35. Crosses show the rotation period data, and short horizontal 
lines the 25th and 90th percentiles of uj, used to characterise the slow and fast rotators respectively. The lines show our models for 
1.0, 0.8, 0.6, 0.42 and 0.28 Mq (respectively), where the solid lines are solid body models, and dashed lines are differentially rotating 
models, with the parameters shown. Plotted are the ONC (1 Myr), NGC 2264 (2 Myr), NGC 2362 (5 Myr), IC 2391, IC 2602 (~ 
30 Myr), a Per (~ 50 Myr), NGC 2516, M34, the Hyades (625 Myr) and the Sun (~ 4.57 Gyr). The IC 2391 data were taken from 
Patten &: Simo nl lll996^ andIC 2602 f romlBarnes et al. (1999 ). The a Per da ta are a compila t ion of t h e results fromlStauffer et al.l 
Stauffer. Hartm ann fc Jonetj ('1989 1). |p7osser et al.l I11993ai') jProssed lll99lD. iPross er et al.' lll993bl) lO'Dell fc CoUier CameronI . 
Prosser et al.l (1995). O'DcU, Hen dry fc Collier CameronI (Il994l). lO'Dell at al.l lll9 96). All ain et al.l (Il996l '). iMarti'n fc Zapatero Oi 



198£) 



1992). 



■ . ■ ■ ■ ■ ■ , , . .. , , - . rici 

1997"), 'Prosser & Randich| (|1998l '). IProsser. Randich fc SimonI l|l998l ). lBarnes et al.l " i(l998i) . and the Hyades data from lRadick et al.l l|l987l ) 
and lProsser et al., (.1995i ). taken from the open cluster database. 



tion (Po ~ 60 days) apparently required at 0.5 Mq is some- 
what implausible in light of measurements of very young 
populations, where typical rotation periods for stars of this 
mass are 10 days. 

In reality, a combination of the two effects is likely, 
which will lead to intermediate values of tdisc and Po. How- 
ever, the required parameter ranges remain implausible in 
light of existing data. It is possible that improved modelling 



of the star u sing a more com plete treatment of differential 
rotation (e.g. ISiUs et al.ir2000l ) will solve this issue. 

6.2.4 Value of LOsnt 

The value of tOsat at each mass is tightly constrained by the 
need to reproduce the NGC 2516, M34 and Hyades data. 
In 16.2.21 we derived values for uisat in coarse mass bins, to 
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Figure 21. As Figure [20l but for varying Pq, with fixed tdisc = 
10 Myr. Tlie curves sliow tlie models for Pq = 1, 2, 5, 10, 20, 
50 days (top to bottom). 



produce model curves providing a reasonable explanation 
for the observations. In order to investigate the mass depen- 
dence of this parameter in detail, we used the NGC 2362 and 
NGC 2516 observations to perform a fit of Wsat as a function 
of mass in much smaller bins to give better sampling, using 
the same percentile based technique. Figure [22] shows the 
results, finding a linear dependence of uisat on mass to be 
a good approximation at the lowest masses, for both solid 
body and differentially rotating models. There is some ev- 
idence that the relation does not apply for higher masses 
(M>0.5 Mq), but in reality Usat is not well-constrained 
here, since the majority of these stars (which are slow ro- 
tators) only reach the saturation velocity for a short time 
during the evolution from 5 to 150 Myr. 

From the fit, at solar mass, 0.8 Mq and 0.5 Mq, the rela- 
tion gives oisat ~ 14.3, 11.1, 6.3 ujq respectively for the solid 
body models, and LUsat ~ 11.9,9.35,5.61 ljq for the differ- 
entially rotating models. The value for solar mass and solid 
body rotation is very close to the value of 14 1^0 we assumed 
earlier, as derived from the convective turnover timescale, 
but at 0.8 and 0.5 Mq, cjsat is larger than the values of 8 ujq 
and 3 luq predicted by Eq. ([TUjl : the slope of Wsat as a func- 
tion of mass is less steep than predicted from the convective 
turnover timescale. 

The decrease of a;sat with decreasing mass, as deter- 
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Figure 22. LOsat plotted as a function of mass, in 0.1 Mq bins, 
from fitting to the NGC 2516 j distribution lower 25 percentile 
(plus symbols) and upper 90 percentile (stars), using NGC 2362 as 
the initial condition, for solid body models (top) and differentially 
rotating models (bottom; Tc = 1.0 Gyr). This diagram indicates 
that a single functional form fails to fit all of these measures 
simultaneously for M>0.5 Mq, but that a simple linear relation 
is a good approximation at lower masses. The straight lines show a 
fit for M < 0.5 Mq, finding {uisat/oJQ) = 16.0 (M/Mq)-1.7 (solid 
body models) and {uisat/(^e) = 12.5 (M/Mq)-0.6 (diff'erentially 
rotating models). 

mined in this work, effectively reduces the amount of angular 
momentum loss experienced by lower mass stars, compared 
to what would be experienced if the value was the same as 
at solar mass. This h as also been seen in previous work, e.g. 
iBarnes fc So"fial lfl99Q ). 

6.2.5 Detailed evolution 

Figure [23l shows the results of attempting to evolve the mea- 
sured rotation periods in NGC 2362 forward in time to the 
age of NGC 2516 using the models we have described, with 
the relations fit from Figure [22] for uisat as a function of 
mass, on an object-by-object basis. By doing this, we can 
test if the model we have presented can reproduce the ob- 
servations of NGC 2516, given the NGC 2362 rotators as 
in input. Note that we have neglected disc locking in this 
analysis for simplicity. 

The apparent mass dependent morphology of the NGC 
2516 distribution is well-reproduced from the NGC 2362 dis- 
tribution by the models, with two exceptions: for the slow 
rotators at the highest masses in the diagram, the predicted 
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Figure 23. Rotation period as a function of mass, using the 
model presented in i|6.2l to evolve the NGC 2362 distribution of 
Figure 1171 forward in time from 5 Myr to 150 Myr for the solid 
body (upper panel) and differentially rotating (middle panel) 
models, and the observed NGC 2516 distribution for comparison 
(lower panel). 

periods are too short by a factor of ~ 2 — 3 compared to 
the observations for the solid body models, and slightly too 
short for the differentially rotating models. This suggests 
that slow rotators experience more angular momentum loss 
on the ZAMS than the models predict. Conversely, the solid 
body models provide a better fit to the fast rotators, espe- 
cially in the mass range 0.6 — 0.8 M0. This is not surprising 
in light of the conclusions of i|6.2.2l 

In particular, the form of mass dependence for 
M < 0.6 M0 is approximately correct, with the exception 
of a small number of objects rotating faster than the ap- 
parent 'limit' of ~ 0.2 days seen in the NGC 2516 data. 



A small number of objects are also visible in the evolved 
NGC 2362 distribution above the P oc relation seen for 
the slowest rotators in NGC 2516. Since these originated as 
relatively slow rotators in NGC 2362, we suspect that they 
are the result of field contamination, or errors in the mea- 
sured rotation periods in that sample due to the limited time 
baseline of the observations. 

This work suggests that the mass dependence of the 
observed rotation period distribution is reproduced by the 
combination of the initial conditions as seen in NGC 2362, 
the mass- dependence of stellar contraction on the PMS, and 
the mass-dependence of the saturation angular velocity ojsat- 
Furthermore, the convergence in rotation rates as a function 
of time (see also Figure [19)) is explained reasonably well by 
saturation of the angular momentum losses, over the mass 
range of interest for NGC 2516, apparently without need to 
invoke any further effects. 

Figure [24] shows the same comparison for M34, showing 
that the feature of constant rotation period as a function of 
mass for Af > 0.6 M© is also well-reproduced by the mod- 
els, at approximately the correct period in the case of the 
differentially rotating models, and again at a shorter period 
than the observed value for the solid body models. It should 
be noted that the adopted values of uosax in this mass range 
are not well-constrained by the NGC 2516 data used to cal- 
ibrate them, and have a strong effect on the morphology of 
this part of the diagram, so it is possible that any deviations 
here result from our assumption of a linear relation over the 
entire mass range. 



7 CONCLUSIONS 

We have reported on results of an i-band photometric survey 
of NGC 2516, covering ~ 1 sq.deg of the cluster. Selection 
of candidate members in &V versus V — I colour- magnitude 
diagram using an empirical fit to the cluster main sequence 
found 1685 candidate members, over a V magnitude range 
of 16 < V < 26 (covering masses from 0.7 M© down to be- 
low the brown dwarf limit). The likely field contamination 
level was estimated using a simulated ca talogue of field ob - 
jects from the Besangon Galactic models (|Robin et al.ll2003l ). 
finding that ~ 650 objects were likely field contaminants, an 
overall contamination level of ~ 39%, implying that there 
are ~ 1000 real cluster members over this mass range in our 
field-of-view. 

From ~ 8 nights of time-series photometry we de- 
rived lightcurves for ~ 104 000 objects in the NGC 2516 
field, achieving a precision of < 1% per data point over 
15 < i < 19. The lightcurves of our candidate cluster mem- 
bers were searched for periodic variability corresponding to 
stellar rotation, giving 362 detections over the mass range 
0.15 < M/Mq < 0.7. 

A striking morphology was found in the rotation period 
distribution as a function of both mass, with a median re- 
lation of P oc (or constant specific angular momentum 
j) for masses below ~ 0.6 M0, with the slowest rotators 
(the majority of the sample) following P oc M'^ {j oc M~^), 
with the fastest rotators over the entire mass range having 
P > 0.25 days (corresponding to a line of j oc M^). Our ear- 
lier M34 results indicate that for M>0.6 M0, the slowest 
rotators follow instead P ~ constant as a function of mass, 
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Figure 24. Rotation period as a function of mass, using the 
model presented in i|6.2l to evolve the NGC 2362 distribution of 
Figure 1171 forward in time from 5 Myr to 200 Myr for the solid 
body (upper panel) and differentially rotating (middle panel) 
models, and the observed M34 distribution for comparison (lower 
panel). 



with a transition between the two regimes at M ~ 0.6 Mq. 
We suggest that this mass corresponds to the transition 
to convectively-dominated stellar interiors moving to lower 
masses. This could be tested by examining rotation period 
distributions over a wider range of ages for M-dwarfs. 

In simple models of the rotational evolution were 
considered, both for the solid body case, and including differ- 
ential rotation between a decoupled radia t ive co re a nd con- 
yectiv e en velope, fo l lowin g iBouvier et al.1 l|l997h and lAUainl 
l| 19981 ) . A iKawaleJ (| 19881 ) angular momentum loss law in- 
corporating saturation at a critical angular velocity tJsat was 
assumed. The results indicate that a mass-dependent satura- 



tion velocity (e.g. of the form suggested bv lSills et al.ll200oh 
is required to reproduce the observations. Furthermore, we 
find that solid body rotation is sufficient to explain the evo- 
lution of the fastest rotators, but core-envelope decoupling 
appears to be needed to explain the evolution of the slow- 
est rotators. This may indicate that the redistribution of 
angular momentum at the core/envelope boundary is more 
efficient in fast rotators. We were unable to find a model 
which fits both the fast and slow rotators simultaneously. 
Further rotation period measurements at low masses, par- 
ticularly for ages intermediate between 5 and 150 Myr and 
for ages > 200 Myr, will provide more stringent constraints 
on the models, but the present observations already indicate 
the need for more theoretical work. 

We intend to publish the final catalogue of NGC 
2516 membership candidates after obtaining follow-up spec- 
troscopy. However, the preliminary catalogue of membership 
candidates is available on request. 
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